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Abstract 

The singularity confinement test is very useful for isolating integrable cases 
of discrete-time dynamical systems, but it does not provide a sufficient criterion 
' for integrability. Quite recently a new property of the bilinear equations appear- 

• ing in discrete soliton theory has been noticed: the iterates of such equations 

o 

are Laurent polynomials in the initial data. A large class of non-integrable map- 
I pings of the plane are presented which both possess this Laurent property and 

have confined singularities. MSC2000 numbers: 11B37, 93C10, 93C55 

3 ■ 

>• . There continues to be a great deal of interest in discrete-time dynamical systems 
that are integrable. There is a vast range of such systems, including symplectic 
maps and Backlund transformations for Hamiltonian systems in classical mechanics 
PP, mappings that preserve plane curves which occur in statistical mechanics, 
discrete analogues of Painleve transcendents 0, partial difference soliton equations 
appearing in numerical analysis and solvable quantum models jl], and equations 

arising in theories of discrete geometry and discrete analytic functions For some 
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time it has been appreciated that integrabihty in the discrete setting is associated 
with certain weak growth phenomena 0, which can be measured by means of suitable 
notions of entropy [7^ or complexity [Hj. 

Given the multitude of application areas in which discrete integrable systems ap- 
pear, the problem of identifying when a given system is integrable is of considerable 
importance. In the continuous setting, the Painleve property has proved to be an 
extremely useful criterion for isolating integrable differential equations 9^, and this 
led Grammaticos, Ramani and Papageorgiou to introduce the singularity confine- 
ment test for discrete equations ^U]. However, while that test has been enormously 
successful at identifying discrete Painleve equations, it turns out that singularity con- 
finement is not sufficient for integrabihty, as was pointed out by Hietarinta and Viallet 
[7]. Those authors found numerous examples of rational maps, taking the form 



which have confined singularities and yet whose orbit structure displays the charac- 
teristics of chaos. This led them to suggest that the stronger requirement of zero 
algebraic entropy (defined in terms of the growth of degrees of iterates) should be 
a necessary property of rational maps that are integrable, in agreement with the 
observations of Veselov [Oj. More recently it was proposed by Ablowitz, Halburd 
and Herbst that the Painleve property can be extended to difference equations using 
Nevanlinna theory [llj, while Roberts and Vivaldi have used the orbit structure of 
rational maps defined over finite fields to detect integrabihty ^2], and Halburd has 
translated the concepts of pJJ into a Diophantine integrabihty criterion for discrete 
equations [T!^ . 

There is a large amount of literature on discrete bilinear equations, including the 
bilinear forms of discrete Painleve equations [H] , and bilinear partial difference equa- 
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tions such as Hirota's difference equation |3j. However, there is one aspect of such 
integrable bihnear equations that researchers on integrable systems have apparently 
overlooked, namely the fact that they have the Laurent property: that is, for suitably 
specified initial data, all of the iterates of these discrete equations are Laurent polyno- 
mials in these data with integer coefficients. It seems that this Laurent phenomenon 
was originally known only to a few people working in algebraic combinatorics, and for 
Hirota's equation it was first proved by Fomin and Zelevinsky within the framework of 
their theory of cluster algebras Pl] , with further combinatorial interpretations being 
found later jT3]. 

One of the simplest examples of a bilinear equation displaying the Laurent phe- 
nomenon was found by Michael Somos, who considered kth order recurrences of the 
form 

[fc/2] 

Tn+k'Tn = ^ ] Tn+k-jT'n+j, k > 4, (2) 

i=i 

taking the initial values tq = ti = . . . = r^-i = 1. Clearly each new iterate Tn+k of © 
is a rational function of the initial data, so one expects r„ G Q, but it was observed 
numerically that for the Somos-4 recurrence 

+ P{rn+2r (3) 
with parameters a = (3 = 1 and starting with four ones, an integer sequence 

1, 1, 1, 1, 2, 3, 7, 23, 59, 314, 1529, 8209, 83313, . . . 

results. Several simple proofs of the integrality of this sequence were subsequently 
obtained (see the article by Gale and the other references ^J). However, it was 
realized that the deeper reason behind this lay in the fact that the recurrence (jHl) has 
the Laurent property, meaning that the iterates are polynomials in the coefficients. 
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the four initial data, and their inverses, and these Laurent polynomials have integer 
coefficients, i.e. r„ G Z[a, l3, t^^, t^^, r^^, r^^] for all n. 

Fomin and Zelevinsky found that their theory of cluster algebras provided a suit- 
able setting within which they could prove the Laurent property for a variety of 
discrete equations ^1], including certain recurrences of the form 



for particular choices of polynomials F (including Somos-/c for k = 4,5, 6, 7), as well 
as integrable two- and three-dimensional bilinear recurrences like Hirota's equation. 
As a parallel development, the connection between the iterates of the general Somos-4 
recurrence (jHl) and sequences of points on elliptic curves has been explained by several 
people: two different approaches are found respectively in the PhD thesis of Swart 
and in the work of van der Poorten ^^1; while independently the author found the 
explicit solution of the initial value problem for both Somos-4 and Somos-5 ^Zj in 
terms of elliptic sigma functions. For earlier unpublished results of Zagier and Elkies 
on the associated elliptic curve and theta function formulae for the original Somos-5 
sequence, see [iH] and [12] respectively. 

An essential observation in the work [17^ was that both the fourth and the fifth 
order Somos recurrences could be understood in terms of a suitable integrable map- 
ping of the plane. For example, setting Xn = r„+ir„_i/(r„)^ transforms Q into the 
rational map 



which preserves the two-form {dxn-i Adxn)/{xn-iXn) and has the conserved quantity 



A symplectic map of the plane with a conserved quantity is the discrete analogue of a 

4 





a (3 



(4) 




Hamiltonian system with one degree of freedom, and hence is integrable in the sense 
of Liouville (see chapter 10 in [20] )■ Thus the map defined by (jH) is integrable and 
belongs to the well studied class |2l E] of rational mappings of the plane that preserve 
an algebraic curve, in this case the curve of genus one defined by (jSj). In fact, certain 
integer sequences satisfying the Somos-4 recurrence ^ were already known to number 
theorists by the name of elliptic divisibility sequences j2l], and these continue to be 
the subject of active research due to the way that new prime divisors appear therein. 

The main result proved below is the following 
Theorem Given a polynomial f{x) of degree d having the form 

fix) = x^' F{x) 

for a non-constant polynomial F with F{0) ^ 0, the recurrence 

Xn+lXn-l = f{Xn) (6) 

possesses the Laurent property if and only if one of the following three cases holds: 
(i) M = and, for all x, the polynomial f satisfies 

fix) = ±x'^ with /(O) = 1 for ^ 2; (7) 

(a) M = 1 and, for all x, the polynomial F (of degree D) satisfies 

F{x) = ±x^ F{l/x), with F(0) = 1 or - 1 for D ^ 1; (8) 

(Hi) M > 2 and F is an arbitrary non- constant polynomial. 

Moreover, up to the freedom to rescale Xn, in cases (i) and (ii) the respective 
conditions (0) and (0) are also necessary and sufficient for the singularities of (0) to 
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be confined immediately, while in case (Hi) the singularity confinement test is failed. 

Clearly to obtain Laurent polynomials it is necessary for f{x) to be a polynomial 
in X, and the recurrence (jH)) is generated by iterating the rational mapping of the 
plane defined by 
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which preserves the two-form 

u = (xy)^^ dx A dy . (10) 

For © to have the Laurent property it is required that Xn & Tl := Z[c, ^, x^^] for 
all n, where c denotes the coefficients of /. To begin with, we take the generic case 
/(O) = A 7^ 0. Given the two initial data Xq, Xi, the next iterates are X2 = f{xi)/xQ 
and X3 = f{x2)/xi, which are both in TZ for any /. The first place where a division 
must occur is for x^ = f{x^)/x2: if X4 G 7^ then /(xi) = X2Xq must divide /(xa). 
Now every Laurent polynomial can be written as a polynomial divided by a monomial 
XqX^, and the ring TZ has the structure of a unique factorization domain, with units 
given by the monomials ±XqX^ for integers i, m. Thus one can do modular arithmetic 
with the elements of TZ in the usual way, and note that xi and X2 are coprime. Then 
X3 = /(0)/xi = \/xi mod a;2 and hence /(X3) = /(A/xi) mod X2, so for divisibility 
of f{xs) by f{xi) we must have that, for arbitrary xi, 

Xf{x,)=^ixif{\/x,) (11) 

where d = degf and the constant fi is the leading coefficient of /. Comparing the 
constant term on each side of (jllll gives = fi^X'^, so /i = iA^^'^/^. If d 7^ 2 then 
to have polynomials in A forces the choice A = 1 (and in any case, the freedom in A 
can be removed by rescahng Xn), so we find /i = ±1 for all d, and thus the necessary 
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condition for the recurrence © to have the Laurent property is that the polynomial 
/ satisfies the reciprocal property ((Tj). 

The condition (jZj) also turns out to be sufficient to ensure that Xn & for all 
n, and also implies that © passes the singularity confinement test. To see that 
this is sufficient for the Laurent property, one can apply Fomin and Zelevinsky's 
Caterpillar Lemma ^Mj , but a direct proof is given here for completeness. Take as the 
inductive hypothesis that Xj ETZfoi < j < n with all adjacent pairs Xj, Xj+i being 
coprime in this range. Then Xn = f{xn-i)/xn-2 = /(0)/a;„_2 = l/xn-2 mod Xn-i, so 
Xn+i^n-i = fi^n) and then using ((Tj) and (jH)) once more gives = /(l/x„_2) = 

±2;;^-2 /(^n-2) = ±x~^2^n-3Xn-i = mod x„_i, whence and Xn+i G 11 

as required. Furthermore, suppose that p is an irreducible element of 71 such that 
p\xn+i] then p\f{xn) and = 1 mod Xn, so p /fx„ and hence x„ and Xn+i are 

coprime, which completes the inductive step. This proves that a;„ G 7?. for all positive 
indices n, and the result extends to negative n by the reversibility of ©. 

As for singularity confinement, note that a singularity can only occur in © if one 
of the iterates becomes zero. So let Xns = a and Xn-2 = r + e where r is any one of the 
(generically distinct) roots of /. Thus Xn-i = /'(?") e + O(e^) ^ as e ^ 0, while 
a;„ = /(0)/r + 0(e) = 1/r + 0(e), which gives Xn+i = (/'(r))-i a e-V(l/r + 0(e)) = 
0(1) because /(r) = implies f{l/r) = by and thereafter all the terms are 
finite (and non-zero) as e — > 0. So we see that the singularity is confined at the first 
possible stage for any mapping of the form (jH)) with /(O) 7^ that has the Laurent 
property. Conversely, it is easy to see that if we start from © with a polynomial / 
such that /(O) 7^ 0, we can always scale so that /(O) = 1, and if we require that a zero 
is confined immediately, then for any root r of /, the reciprocal value 1/r must also 
be a root, which implies that (|7j) must hold, and hence the mapping has the Laurent 
property. 

In the above considerations we imposed the restriction /(O) 7^ 0. We now describe 
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the situation for the case /(O) = 0, and so set f{x) = F{x) with integer M > 1, 
-F(O) 7^ and degF = D > 1. There are two further cases to consider, according to 
whether M = 1 or M > 2. In case (ii), we find that for the Laurent property to hold 
we must have 

f{x) = X F{x), F{x) = ±x^ F{l/x) (12) 

for all with -F(O) = 1 or — 1 for D ^ 1. So F must satisfy the same reciprocal 
property as for / in the generic case, and once again this condition is also sufficient; 
the proof of the Laurent property in this case is slightly more involved and will be 
presented elsewhere. Similarly to the argument for (0), it is easy to verify that the 
condition (fT^ also implies that the mapping (jHl) passes the singularity confinement 
test. The other case (iii) is somewhat different, for upon taking 

/(x) = M>2, F arbitrary, (13) 

we find that the mapping defined by (jH)) gives x„ G 7?. for all n. To see this, take 
as the inductive hypothesis that x„_i, x„ and p„ = (x„/x„_i) G TZ, and then write 

Pn+l = (Xn+l/Xn) = Pn X^^^"^ f {Xn) E TZ hj the liypothcsis, so Xn+1 = Pn+l^n € TZ 

as required. However, to see if a zero is confined in such a mapping we set Xn = a, 
Xn+i = r + € where F(r) = 0, so that Xn+2 = 6*26 + O(e^) ^ as e — > 0, where 
C2 = a^^r^' F'{r). Subsequent terms have Xn+z ~ C^e^'', Xn+4 ~ Cit^^^^^, Xn+5 ~ 
Cse*''^^"^*^ etc. for certain (non-zero) constants Cj, so the powers of e continue to 
grow and the zero is not confined. (However, note that we have explicitly excluded the 
trivial case F = (3 = constant, f{x) = jSx'^'^ , which always has the Laurent property 
and satisfies singularity confinement in the sense that a; = cannot be reached from 
any non-zero initial data.) 

We have seen that there is a close connection between the Laurent property and 
singularity confinement for discrete equations of the form (jH}, but what about the 
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integrability of such mappings? For a measure-preserving mapping of the plane to be 
integrable it must have a conserved quantity. We can start by considering the lowest 
degree examples of polynomials /. If ci = 1 then 



Xn+lXn-l = 'JXn + S, (14) 

and any such map is integrable because it has a conserved quantity K defining an 
elliptic curve, i.e. K equals 

Xn-l+Xn^ . (15) 

Xyi—lXri 

Equation (fT^ is bilinearized via x„ = / ijn+i'Tn) , yielding the special Somos-7 

recurrence 

which can be related to a Somos-5 recurrence and thence solved in terms of elliptic 
sigma functions [TTl . However, if we require that ()14|1 itself should possess the Laurent 
property, then either 5 = with 7 arbitrary and the map cycles with period 6, or we 
are in the situation (0) so that 7 = ±1 with 5 = 1 and the map cycles with period 5. 
It is interesting to note that in the latter case (fixing 7 = 1) this map is equivalent 
to the functional relation that appears in the thermodynamic Bethe ansatz for an A2 
scattering theory [22] • 

The case d = 2 is special: if A = /(O) 7^ then the Laurent property requires ((Tj) 
to hold, and there is the extra freedom to leave A 7^ arbitrary. Then there are two 
choices of sign, giving either 

a;„+ix„_i = + i/x„ + A, (16) 
or alternatively Xn+i Xn-i = —x"^ + A. The second choice can be transformed into the 
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first by taking x„ = /€„ x„ where = 1 and = — 1 for all n, so we can just 

consider p6|) which has a conserved quantity that defines a conic, namely 

Furthermore, the iterates also satisfy a linear recurrence of the form (Q), viz. 

Xn+l + Xn-l = LXn — 

SO in terms of the initial data we have L = L{xo, Xi) G 71 and also Xn G Z[z/, xq, Xi, L] C 
7?. for all n, which is even stronger than the Laurent property. When d = 2 and 
/(O) = there are further sub-cases. If (fT^ holds then (up to rescaling x„ u'^Xn) 
either f{x) = x{x + l), which is just a special case of (fTB|) . or we have the opposite 
choice of sign and the mapping is given by 

Xn+l Xn-l = -Xn{Xn + 1), (18) 

which is also integrable, having the conserved quantity 

J -^n— 1 ~^ ■^n ~^ {Xn—1 (2x„_i + 2x„ + 1) 

~ 2 9 ■ 

■^n-l-^n 

Moreover, for (|18|) the iterates also satisfy the sixth order linear recurrence 

Xn+Q -\- {J — l)(a;n+4 — Xn+2) — a;„ = 0, 

which provides an alternative proof that Xn & based on the the fact that J E TZ 
and Xj G 7?. for j = 0, 1, . . . , 5. 

For (i > 3 all of the maps © with the Laurent property are non-integrable. To see 
this, one can count the growth of degrees to show that the algebraic entropy is non- 
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zero 0, but a simpler test to apply is Halburd's Diopliantine integrability criterion 
|13j . For a rational number x = p/q the logarithmic height is h{x) = logmax{|p|, |g|}, 
and Halburd's test requires that for rational- valued maps to be integrable, h{xn) must 
grow no faster than a polynomial in n. Suppose that a sequence of (real or complex) 
iterates is such that — oo and A„ = log ~ C^" for real C > 1 and some C > 0. 
Then taking logarithms of both sides of © gives A„_|_i + A„_i — dAn ^ and hence 
with (i > 3 we find ( = {d+\/d? — 4)/2 > 1 as required. Now for equations © having 
the Laurent property, if we take / to have integer coefficients and set Xq — X\ — 1 
then all x„ are integers, so that A„ = h{xn)- Moreover, if all the coefficients of / are 
positive then the terms of the integer sequence will have precisely these asymptotics, 
so that the logarithmic height grows exponentially and lim„_>oo(log/i(x„))/n = log( 
(which also happens to be the value of the algebraic entropy for these maps). 

To understand the deep connection between singularity confinement and the Lau- 
rent property, we propose to extend the above results in at least three directions. 
Firstly, given a pair of polynomials /i, /2 each satisfying ((Tj), it is simple to prove 
that the composition (pi ■ 02 of the corresponding pair of maps also has the Laurent 
property; the choice fi{x) = x^ + 1, f2{x) = x'^ + 1 generates a cluster algebra of 
rank 2 [SB] • Secondly, there are many higher order discrete equations (integrable and 
non-integrable) with the Laurent property. Thirdly, this property can also apply to 
non-autonomous equations, such as the bilinear forms of discrete Painleve equations. 
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